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■ Abstract. Caustics are formally singular structures, with infinite den- 
sity, that form in collisionless media. The non-negligible velocity dis- 

■ persion of dark matter particles renders their density finite. We eval- 
^ ' uate the maximum density of the caustics within the framework of 

VO , secondary infall model of formation of dark matter haloes. The result 

is then used to demonstrate that caustics can be probed by properly 
' stacking the weak-lensing signal of about 600 haloes. CFHTLS accom- 

panied by X-ray observations and the space-based experiments like 
SNAP or DUNE can provide us with the required statistics and hence 
a way of distinguishing between the viable dark matter particle candi- 
. dates. The extension of our results to more realistic models including 

' the effects of mergers of haloes is briefly outlined. 

O 



o 

in 
o 



1 Introduction 



Under the gravitational instability, cold streams of dark matter with different ve- 
. ^ \ locities cross at caustic surfaces. The density of dark matter at caustics formally 
^ diverges and in the regions bounded by caustics the velocity field is multi-valued. 
^ The full description of the evolution of the fluid is provided by the Vlassov-Poisson 
equation, in 6-dimensional phase space. Direct numerical integration of this equa- 
tion remains a challenging task and so far is best achieved for systems with low 
number of dimensions and or with strong symmetry constraints (see e.g. Alard 
& Colombi 2004 and references therein). The most popular alternative to the full 
integration of this equation is to sample the phase space with discrete particles. 
However, extremely large number of particles is needed to sample the fine- structure 
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of the phase space and significant softening of the forces is required to avoid spu- 
rious collisional effects which would smear out the caustics (see e.g. Binney 2004, 
Melott et al. 1997). 

Until a full 3-dimensional solution of Vlassov-Poisson equation or ground- 
breaking resolution in N-body simulations or other innovative solutions are achieved, 
(semi-) analytic modeling of haloes can provide approximative description of the 
dynamics at fine scales. 




Fig. 1. Left panel: Phase space diagram (radial velocity versus radius) of a Milky- Way- 
type halo, taken from a PM simulation with 16 millions of particles on a 1024 grid in a 
4 Mpc simulation box realized at IDRIS computing center (CNRS). The halo has about 
8 million particles in its virial radius and contains only 40 substructures, as they are 
rapidly disrupted due to softening of the forces in the PM routine. The smoothness 
of the gravitational potential (there are many particles per softening length) enables us 
to see the caustics since the phase-space structure can survive the relaxation effects. 
Right panel: Secondary infall model's prediction of the phase space of a spherical halo in 
smooth accretion. The caustics in the model are far better resolved as compared to the 
left panel, due to the simple assumption of sphericity, no abrupt major merger and no 
substructures. The panel clearly shows that the streams cool down as they collapse: a 
consequence of phase-space volume conservation as filamentation builds up. 



Analytic works on the halo density profile, and its caustics, started with the 
works of Gott (1975) and Gunn (1977). In an Einstein-de Sitter Universe a spher- 
ical over-density expands and then turns around to collapse. After collapse and 
at late times, the fiuid motion becomes selfsimilar: its form remains unchanged 
when length are re-scaled in terms of the radius of the shell which is currently 
turning around (rta) and falling onto the galaxy. Self-similar solutions give power- 
law density profiles which is convolved with many small-scale spikes i.e. caustics 
(Fillmore and Goldreich 1984, Bertschinger 1985a, 1985b). 
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In the self-similar accretion model, Newton's law is reformulated as 
£X 7dX 8 ^ 2 

by using the non-dimensional variables 

A=4\ ; C = ln(±] ; M(A)- 



where m(r, t) is the mass inside a radius r [m(r, is not constant due to shell- 
crossing], ph = l/GnGt'^ and tta is the turnaround time for a given particle (i.e. 
when the particle is at its largest radius). This reformulation also assumes the 
power-law initial condition 

SMi ( Mi 



- 



(1.2) 



where Mq is some reference mass and Mi is the unperturbed mass within the initial 
radius , It might appear that the asymptotic self- similarity is a direct outcome of 
the power-law initial mass profile. However, it has been shown that starting from 
a generic initial condition, a power-law density profile is achieved prior to the first 
turnaround time (Moutarde et al 1991). The turnaround radius increases with time 
as rta ^ t'^2/3)(i+i/3e) hcucc the mass inside the turnaround radius increases 
as Mta ~ (1 + 2;)"^/^. This scaling compares with that of the characteristic mass 
in a scale- free hierarchical clustering Universe ?^ (1 + 2^)"^/^^^+^^ where the 
power-spectrum of initial density fluctuations is P{k) ^ k^. A central point mass 
perturbation corresponds to e = l,rta ~ t^/^;Mta ~ 1/(1 + z) which is the only 
case considered here. Extension of our results to more realistic initial conditions 
{e.g. e ^ 0.2) will follow in forthcoming works. 

For a perfectly cold dark matter medium, the density profile close to the kth 
caustic at is (Bertschinger 85b) 

(The values of the various quantities ^/c, A/c, A^' are given in Bertschinger 1985b.) 

When the temperature of particles is not strictly zero, caustic positions are 
shifted by a small value ^A and the density near the caustic is modified as: 
Per (A) = J dv po[X — SX{v)] f{v) . A simple top-hat velocity distribution function, 
f{v), yields 

P.(A) ^( TT^ 
PH V A,A^ i 



K-^ -^K->^ for A<A^, 

A^ - A for A~ < A < A^ , 



for A > A+ , 



(1.4) 
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for the density near the k-th caustic, where = — Aj^ and = + Aj. 
and the thickness of the k-th caustic is A/^. A simple re-scaling demonstrates that 
these densities are universal: they are valid for all caustics and irrespective of size 
and mass of their halo, (for full details see Mohayaee & Shandarin 2005). 

Caustics could have immediate impact for dark matter search experiments. 
Their high density and the fact that they stay well-seperated from each other (see 
left panel of Fig [2j can leave detectable fluctuations in the annihilation products 
of the dark matter particles [if it constitutes for example of self-annihilating axions 
or neutralinos] (see e.g. Sikivie & Ipser 1992, Sikivie et al. 1997) and or in the 
weak lensing data, which we shall study here [see i.e. Hogan 2001, Charmousis 
2003, Gavazzi et al. 2005, Onemh 2004 and 2005]. 

For weak lensing, the Abel integral relates the 3-dimensional density (p) and 
the 2-dimensional density profiles (S) by (see Gavazzi et al. 2005 and references 
therein for full details) 

We numerically integrate the above expression, for the density profiles (|1.3|) and 
(|1.4|) . We define a pseudo-shear: r(A) = (S — T?)/ pnrta and the corresponding 
noise level Vn- For an EdS cosmology, and considering an annulus of inner and 
outer radii Ai and A2 respectively, it is straightforward to write Fat in the following 
units: 



_ ^ X Dos /5Mpc\^, ,0 /30arcmin f (Je\ 1 

r„(A.,A.) = 2.10- (^) (1 + ,.f X ^ __ (-) . 

(1.6) 

Fig. [2I shows pseudo-shear as a function of distance from the center of the clus- 
ter for the thermal velocity dispersions a = km/s and a = 60 km/s. Comparing 
these curves, one can see that the sawtooth patterns due to caustics survive sig- 
nificantly high temperatures. Next, we consider the noise level for a fiducial halo 
at redshift zi — 0.3 [since intermediate redshift haloes at zi ~ 0.2 — 0.5 are the 
most likely targets; see Gavazzi et al. 2005 for further details], and a turnaround 
radius rta = 5Mpc which is a typical value for clusters (upper green binned curve). 
With a single halo the detection of caustics is impossible. Since, the noise falls 
with the square-root of the number of lenses (as well as the source), if one could 
stack the signal from a few tens of clusters, the noise level will be low enough to 
be sensitive to caustics as a whole (lower blue curve). Systematic error can arise 
due to the intrinsic ellipticity of the clusters. The precise effect of which can be 
well- studied with the help of N-body simulations in future works. However, due 
to its significant separation from the second caustic and the turnaround radius, it 
is likely that the first outer caustic would survive smearing due to substructures 
and ellipticity of the haloes. 

Although in this work we have considered dynamics of a simple isolated haloes, 
we believe that the results can be extended to haloes growing in cosmological en- 
vironment. One of the short-comings of the secondary infall model could be that 
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Fig. 2. Left panel Left panel ; Main plot The top line (filled circles) shows the separation 
of caustics, Afc — A^+i, as a function of the distance, Afc, from the center. The bottom 
line shows the thickness of the caustics. Left panel ; Inset: The ratio of the thicknesses 
of the caustics to their separations , 26Xk/{Xk — Afc+i), is shown as a function of their 
radii. This panel demonstrates that in the course of the gravitational evolution streams 
remain well-isolated from each other in spite of the fact that their separations diminish. 
Right Panel: Pseudo-shear is plotted for two values of a: cold medium a = (black 
curve) and warm medium a 60 km/s (red curve). Unphysically large value of velocity 
dispersion is therefore required to smear out the caustics. The green (resp. blue) binned 
curve is the noise level for one (resp. 100 stacked) halo(es). 




Fig. 3. Phase space of a ID- ID simulation of merger of two haloes using the waterbag 
code of Colombi and Touma (2005). The simulation indicates that in the merger of two 
haloes, phase-space folds do not mix and the topology of caustics remains intact. 



it ignore the hierarchical scenario of merger. Merger might lead to distorted topol- 
ogy for the caustics and hence render them practically non-detectable. However, 
minor mergers which are on average isotropic are well-described by the self-similar 
model, which assumes a spherically symmetric accretion. Thus such mergers would 
not erase the caustics. Indeed even if the condition of isotropy is not satisfied, the 
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caustics retain tiieir siiapes quite well. The numerical solution to a 1-dimensional 
Vlassov-Poission equation clearly suggests this fact in Fig. [31 In mergers of un- 
equal mass haloes the smaller halo wraps around the larger halo without disrupting 
its phase-space structure. In fact it populates it in a coherent manner as shown 
in the top panel of Fig. [HI 

In a major merger, a similar situation follows. The phase sheets of the two 
caustics do not mix, However, they fold closely together in phase space [bottom 
row of Fig. [SJ. Thus, Fig [31 demonstrates that in the merger of two haloes the 
topology of caustics remains intact. 

Another important feature demonstrated well in Fig. |31 is the gravitational 
cooling effect. With the increase of dynamical time (t) the coarse-grain velocity 
dispersion which is related to the number of folds in the phase space increases. 
However, each fold cools down and becomes thiner, a surprising feature of a sys- 
tem with negative specific heat, but which can be understood roughly by the 
conservation of phase-space volume. 

For a perfectly cold dark matter with zero velocity dispersion, caustics are 
singular objects with zero thickness and infinite density which form inevitably in 
the course of gravitational collapse. The small velocity dispersion of dark matter 
particles renders the density and the thickness of the caustics finite. The latter 
quantities, if observed, could provide direct evidence for the existence of dark 
matter and put bounds on its mass. 

We have shown that the existence of dark matter caustics could be probed by 
properly stacking the weak lensing signal of a reasonable number of haloes. The 
main observational limitation is perhaps the precise estimation of the turnaround 
radius, rta, of superimposed haloes. However, the loss of a few per cents relative 
accuracy in the determination of rta (or asphericity) can be compensated for by 
stacking about 600 haloes. Wide field surveys such as the ongoing CFHTLS accom- 
panied by X-ray observations can provide the required statistics for a successful 
detection of caustics. The number of haloes required to be superimposed will be 
lowered by a further factor of 3 for future space-based experiments like SNAP or 
DUNE. 
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